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Abstract: Asymptotic methods in nonlinear dynamics are used usually to improve 



perturbation theory results in the oscillations regime. However, for some problems 
of nonlinear dynamics, particularly in the case of Higgs (Duffing) equation and the 
Friedmann cosmological equations, not only small oscillations regime is of interest but 
also the regime of rolling (claiming), more precisely the rolling from a top (claiming 
to the top). In the Friedman cosmology, where the slow rolling regime is often used, 
the rolling from a top (not necessary slow) is of interest too. 

In the present work a method for approximate solution to the Higgs equation in 
the rolling regime is presented. It is shown that in order to improve perturbation 
theory in the rolling regime it turns out to be effective not to use an expansion in 
trigonometric functions as it is done in case of small oscillations but use expansions 
in hyperbolic functions instead. This regime is investigated using the representation 
of the solution in terms of elliptic functions. An accuracy of the corresponding 
approximation is estimated. 
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1. Introduction 



Known asymptotic methods in nonlinear dynamics such as the Krylov-Bogoliubov 
averaging method, the Lyapunov method of slow varying parameters, the Poincare 
and the Van der Pol methods, and KAM-theory are used to improve perturbation 
theory results when describing small oscillations, see [J| |[ |[ |J . Even the term itself - 
nonlinear mechanics - is frequently considered as a synonym to nonlinear oscillations 
theory, see 0. 

However, in case of the Higgs equation in field theory [EJ, the Friedman equation 
in cosmology [§], 0, § and some other problems of nonlinear dynamics not only regime 
of small oscillations is of interest but also the rolling regime, see for example P, IT] 
and references within. 

Small perturbations of quasi-periodical motions usually are considered as per- 
turbations to the linear equation of the harmonic oscillator: 

q + co 2 q = 0, (1.1) 

where q = q(t) is a real- valued function of time t and u > 0. 

Rolling (or climbing) regime by definition is a small perturbation of the following 
linear equation: 

g-A = 0, (1.2) 

where /i > 0. 

In the present work a method to study nonlinear systems in the rolling regime is 
presented. This method can be considered as an hyperbolic analogue of the averaging 
method 0. We consider it in an example of the Higgs equation 

q\t)-ii 2 q{t) = -eq 3 {t), fi > 0, e > 0, (1.3) 

where e is small coupling constant. The method gives approximate solutions using 
expansions of the exact solution in hyperbolic functions. The following approximate 
solution is found in Theorem 1 in Sect. 7.3: 

/ / 3eA 2 \\ 1 e 2 A 3 
qapprit) = Asinhf tfil 1 + -— J J - — -^sinh(3t/i) , (1.4) 

where A is an arbitrary constant. One has the bound 

\q(t) - q appr (t)\ < Ce 2 -° , < a < 3/4 (1.5) 

which is valid for time in the interval 

0<t<cilog— (1.6) 
e 
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that expanding with the vanishing e. 

It is well-known that simple equations of nonlinear dynamics allow exact solu- 
tions in terms of elliptic functions. Correspondence between these exact solutions 
and asymptotic averaging method is considered in the book 0, where an equation 
for anharmonic oscillator with real frequency is investigated 

q(t) + u 2 q(t) = -eg 3 (t), u > 0, e > . (1.7) 

In fH numerical comparison of approximated and exact solutions is presented. 

In the present work, unlike ||, the equation ( |1.3|) with the imaginary frequency 
is considered. It also has a representation in terms of elliptic functions. The rolling 
regime for solutions to this equations is considered. Equations ( |1.7|) and ( |L3l) are 
connected via change u = ifi. However, the averaging method is not applicable to 
study the rolling regime for solutions to ( |1.3| ). We use another method that con- 
sists in expansion of solutions in hyperbolic functions. 1 Applications of approximate 
solutions in functional mechanics are considered in ]TTl [T2], O . 

It is well known that by a shift at a constant q(t) = x(t) + n/y/e, the equation 
( P-.3Q is converted into the equation describing nonlinear oscillations 

x + 2/j, 2 x = — ex 3 — 3y/efxx 2 . (1.8) 

However this transformation does not allow to investigate the rolling mode more 
effectively then explicit consideration of the initial Higgs equation. 

The paper is organized as follows. In Sect. 2 we discuss a hyperbolic analogue 
of the averaging method. In Sect. 3 solutions to the Higgs equation in terms of 
elliptic functions is presented. Expansions of elliptic functions in hyperbolic ones are 
presented in Sect. 4 where bounds on the n-mode approximations are also presented. 
A rearrangement of the expansion of the elliptic function in the form suitable to make 
connection with the series obtained via the hyperbolic analogue of the averaging 
method is exhibited in Sect. 5. The identity of these two series is noted in Sect. 5. 
Estimations of precision of hyperbolic approximations are presented in Sect. 6. In 
Appendix some details of the proofs of Lemma 1 and Theorem 2 are presented. 

The paper is based on |T7|, [LB]. 



: It would be interesting to consider possible correspondence of the method of expansion in hyper- 
bolic functions used in the present work with small t with the first Lyapunov method development 
Q that is effective with t —} oo. 
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2. Hyperbolic analogue of the averaging method 

The simplest equation considered in nonlinear dynamics has the following form: 



q + u 2 q = ef(q,q), (2.1) 

where q = q{t) is a real- valued function of time t , frequency a; is a positive number, 
e is a small parameter. 

In applications the following equation is also considered 

q-fi 2 q = ef(q,q), (2.2) 

where /i 2 > 0. It comes from (2~I) via a change u = i\x. This equation is called an 



equation of rolling. For instance, the Higgs equation has this form (see 0]) and also 
the Friedman equation with a tachyon field in cosmology (see [IIJ ) . 

By analogy with the Krylov-Bogoliubov expansion method, used to solve equa- 
tions describing oscillations fl2.1|) , we will look for general solution to the equation 
describing rolling ( p. 2D in the form of expansion 



q = asmhip + eui(a,ip) + e 2 U2(a,ip) + ... (2.3) 

Here a, ip as functions of time that are determined by differential equations 

a = eA(a,e) = eA l (a)+e 2 A 2 (a)+...,tp = fi+eB(a,e) = n+eB 1 (a)+e 2 B 2 (a)+... (2.4) 

Let us consider equation 

q(t) -fi 2 q(t) = -eq 3 {t), fi > 0, e > 0. (2.5) 

We look to solution of the form 

q(t) = 21 sinh(i/im) + ...., (2.6) 

and define parameters 21, m in a way to keep the functional form of the solution after 
substitution in the equation (|2.5|) . We get on the right-hand side of ( p.5| ): 

- e(2lsinh(t / um) + ....) 3 = -e ^2l 3 sinh(3t//m) + 3 sinh(t//m)) + ...J (2.7) 

On the left-hand side of ( |2.5| ) we obtain: 

{d 2 - ^ 2 )(2tsinh(t//m) + ....) = 2l/i(m 2 - 1) sinh^m) + .... (2.8) 
We see that in order to get an equality we have to equate 

2l/i 2 (l-m 2 ) = -e^2l 3 , (2.9) 
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Therefore it follows that 

9 3e2l 2 , 

- 2 = i + i-r- (2-io) 

In order to keep the form of the approximation in the right-hand side and the 
left-hand sides of ( |2.5| ) one has to take the next term of the approximation in the 
form 

3e2F 



sinh ^t/i^l 



q(t) = 2lsinh( tfi[l + -—^ ) ) + Sti sinh(3t/i) + ... . (2.11) 



Taking into account ( |2.7|) we get 



1 e 2 2l 3 

2li = - — ^4-. (2.12) 
32 /i 4 v ; 

Approximation ( p. ll|) ,( pTT2|) is considered in ||. 

The described procedure can be carried on and it turns out that for equation (|2.5|) 
the procedure leads to the representation of the solution with initial data q(0) — 
in the form of an infinite sum of n-mode terms 

oo 

?(*) = (2-13) 



n=0 



where 



s (t) = 2lsinh(//mt), m = m(A), (2.14) 

*„(*) = 2la n (A)sinh((2n + l)//mi) , n = l,... (2.15) 

Here m (A) and a n (A) are power series of dimensionless parameter 

2l 2 e 

A = (2.16) 



m(A) = l + miA + ... +m„A n + ... 

a„ (A) = a n0 A n + a„i A n+1 + ... + a„„ A 2 ™ + ... 



(2.17) 
(2.18) 



The solution depends on parameters e, fi, and it is specified by the only parameter 
21. Within this notation the first terms of fl2.13|) are written in the form 



q(t) = 2lsinh ( ji[ 1 + 
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e2l 2 



8/x 2 
21e 2 2l 4 



1521 4 

2V 



(2.19) 



32/x 2 102V 



sinh( 3/x( 1 + ^e It 



21e 2 2l 5 
102 V 



sinh(5/xi) + 



This series can be formally obtained for the result of the averaging method (see (|3.7Q 
below) by the change u — > ifi, A — > — i2l. 
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For the Higgs equation with friction: 

(<9 2 + 2ehd - tu 2 )q + eg 3 = 0, q(0) = (2.20) 
an approximate solution has the form 

q(t) = 2te-^sinh Lft - ^-(1 - e~ 2 ^))) (2.21) 
2l 3 _ 3eht 32l 2 (l-e- 2 ^^ 



- e 3V e Sinh ^^ Toy^JJ ^ 

The method developed in this work can be also used to investigate a rolling mode 
(not necessarily slow one) for the Friedman equation either in local or in nonlocal 
cosmology, see M, Pfl: 



/ '2 2 

q + e lS j |- ± a; 2 |- + e 2 y(g) ± o; 2 g + e 2 K'(g) = 0. (2.23) 

3. Higgs equation and elliptic functions 
3.1 Solution to the Duffing equation 

In this section the known solutions to some nonlinear equations are presented. 
As it is known the equation of the form 

q + V'{q) = 

has the integral of motion: 

E = \q 2 + V(q) 

and solution to the equation can be presented in the form of quadrature 

dq 



t 



'2{E-V{q)) 



If V(q) is a quartic polynomial, then the last integral can be represented in terms 
of elliptic functions @, pf, [||. 



It follows that the solution to the equation for the anharmonic oscillator (the 
Duffing equation) 

q{t) + u 2 q{t) = -eg 3 (t), u > 0, e > (3.1) 

has the representation 

q(t) = acn(nt + b,k), (3.2) 
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where cn(u, k) is the elliptic cosine function of argument u and modulus k, a is the 
amplitude and b is the phase. The frequency Q and the modulus of elliptic function 
k are obtained from parameters of the equation u, e and depend on amplitude a: 



«=v^W, * = ^£-5j. (3-3) 
The function cn(Qt + b, k) can be expanded in trigonometric series 

27T A Cf- 1 ^ 



n=l x 7 



(3.4) 



where 

q = e -w T^, q'^e - ^. (3.5) 

Here K = K(/c) is the complete elliptic integral of the first kind and K' = K(fc'), 
where fc 2 + A;' 2 = 1 [ 14| , ||15|| , [|T~6]| . In particular, with 6 = — K one obtains: 



„ oo re— 1/2 / 0+\ 

cn(« - K, *) = £ £ ^ sin((2 n - I,**) (3.6) 
From this representation one can obtain an asymptotic expansion for small e: 



here A is an arbitrary real parameter. This expansion one can also obtain by means 
of the Krylov-Bogoliubov averaging method 0. 

3.2 Solution to the Higgs equation 

Solution to the equation 

g- /i 2 g=- eg 3 , e > 0, (3.8) 

can belong to one of the three types depending on initial conditions 

9(0) = go, q{0)=v o . (3.9) 
A type of the solution is defined by energy 

E = \vl - l -n 2 ql + \eql (3.10) 
Solutions have one of the following forms 
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Figure 1: Potential energy and initial conditions. . The initial condition corresponding to 
a movement in two holes, E > 0; B. The initial condition, corresponding to a movement in 
two holes with infinite period, E = 0; C. The initial condition corresponding to a movement 
in one hole, E < 0. 

• E > 0; the two-holes solutions (Fig 1A) 

q(t) = acn{Qt + b,k) (3.11) 



jjt I I AeE 



« 2 = -(l + A/l + ^f 1 (3.12) 



^ 2 = /A/l + ^ (3-13) 

v + ~ 

. 2 1 1 1 



2 2 /j + 46E 



(3.15) 



Parameter b is defined from the condition 

q(0) = acn(b,k). (3.16) 

We will discuss the case with initial data of the special form g(0) = 0. In this 
case 

b=-K (3.17) 



and 



For this solution 



q(t) = acn(nt - K, k). 



VI - 2k' 2 

The energy E is expressed in terms of parameter k', 
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yU 4 "~ (2 k' 2 -I) 2 

We also note the dependency between a, e and kl 



e V 1 -2k' 2 

and Q with \i and k! 



a 2 = ^(l + —l—\ (3.19) 



Q 2 = p - - (3.20) 
1 - 2k' 2 K ' 



E = ; the two-holes solution with period T = oo 



»<«>="Vis5(ibi) (3 - 21) 

Parameter b is defined by initial condition 

V e cosh(o) 

< 0; the one-hole solution, for instance, the solution in the right one 

q(t) = adn{Qt + b,k) (3.23) 



fc 2 = 2-^ = 2- 7 (3.26) 

(i + yrr^) 

Parameter b is defined by the relation 

q(0) = adn(b,k) (3.27) 

In the case with initial data of the special form v o = we have 

q(t) = adn(m-K,Jfe). (3.28) 

The rolling mode corresponds to a movement in the vicinity of the turning 
point, that is closest to the top of the potential (see Fig.[L] C). 



- 9 - 



4. Elliptic functions expansion 



We consider solution ( |3.18 ). This solution is periodical with period T = 2K/1X Here 
K = K(A;). We show that for t, such that \t\ < T/2, the solution can be expanded 
in hyperbolic series (|2.19|) . We also show that the first terms of the expansion ( j2.13|) 
do represent an approximation to the exact solution ( |3.11| ) with b = — K. 

This remark is of importance because in a case of more complicated potentials 
or in a case of a problem with the friction we can explicitly apply the hyperbolic ana- 
logue of the Krylov-Bogoliubov method, while the exact solutions remain unknown, 



see, for instance, the problem with friction, (|2.20 ). 

We need an expansion of elliptic functions about the point that is equal to minus 
half-period, i.e. cn(u — K, k), for small k'. In this case K — > oo as 



v , 4 k '\ 4 k ' 2 
K ^ In — H In— h 

k' 4 k' 2 4 



According to (|3.15 ) if e — >■ 0, then k' — > and vice versa, 



(4.1) 



k' 



The summation formula for the elliptic functions together with 



(4.2) 



sn(K) = 1, cn(K) = 0, dn(K) = k' 



(4.3) 



gives 



cn(u — K) 



cn(w)cn(K) + sn(ti)sn(K)dn(-u)dn(K) 

1 — &; 2 sn 2 (w)sn 2 (K) 
sn(u)dia(u)k' sn(u)dn(u)k' ,sn(w) 

= th 



A; 2 sn 2 (w) 



dn 2 (M) 



dn(w) 



(4.4) 



4.1 Perturbation theory 

In this section we investigate an expansion of the solution to (|2.5|) in terms of the 
elliptic function in the case of small parameter k! (i.e. e) and compare the first terms 
of the series with the exact solution. 



For small k' the following expansions take place ||15|| , p. 386, equations 16.15.2- 
16.15-3, 



sn(w, k) = tanhu H — k' (sinhwcosh-u — u) 5 1" ••■ 

4 cosh u 



dn(w, k) 



1 



1 



cosh u 4 
Consequently, for small k' we have 



-k (sinh u cosh u + u) 



sinhw 
cosh 2 u 



+ .. 



(4.5) 
(4.6) 
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Sn(u,k) 1 2 . . . 2 v 

-— — - = smn« + -K smhu- sinh w) 
dn(u, fc) 4 v ; 

l k' 2 u( — \ tanhu]+... (4.7) 



4 \ cosh it 
and we get the following expansion 

cn(w - K, k) = k'c {u) + F 3 Ci(m) + ... (4.8) 

where 

Co (it) = sinhw (4.9) 

c-i(u) = — (sinh u — sinh 2 u) u ( — : tanhw | (4.10) 

v ' 4 v ' 4 V cosnM / 

Expansion ( |4.10| ) contains a secular term k' 3 u, however, unlike trigonometric 
case, it is dominated by the term ^k' 3 sinh 2 u with large u. First two terms of 
the expansion ( 4.10|) are presented in Fig.0 A. Modified perturbation theory will 



be developed using an expansion in hyperbolic functions in the next subsection. 
This modified theory gives a better approximation to the exact solution for large 
u. It is interesting to mention, that, as it will be shown in the next subsection in 
the modified perturbation theory the term — ^A: /3 sinh3u contributes greately. Its 
contribution increases faster with increase of u, than the term from perturbation 
theory jk' 3 sinh 2 u. 

In the Fig. || the comparison of the first two terms of perturbation theory with 
exact solution is presented. The first approximation is the upper solid line, the sec- 
ond one is plotted by the dotted line. The lower thick line shows the exact solution. 
It can be seen that the first and the second approximations of conventional pertur- 
bation theory almost coincide with each other. The second approximation calculated 
according to the conventional perturbation theory is not in fact an improvement of 
the first one. 

4.2 Modified perturbation theory 

We will show that the function cn(w — K, k), with \u\ < K and small k', allows the 
following representation 

oo 

cn{u — K,k) = ^T(-l) n An sinh((2n + 2)u) (4.11) 



where 



a ...... . D 

Figure 2: The first two approximations comparison, calculated according to conventional 
perturbation theory, (4.9) and (4.10), with the first two approximations of the modified 
perturbation theory, (4.15) and (4.16), and with the exact expression cn(u — K, fc), k = 
0.999. 



A 



71 



1 



kK' cosh)(2n + l)p') 



(4.13) 



and the series (|4.11 ) converges with \*JKu\ < K and 

In order to prove ( 4.11|) we use 
tation holds, see equation (133), p. 85 



-l< k < l. 

. Note that for the following represen- 



sn [u) 



7T 



J sinh u' sinh 3u' sinh bp' 
dn(w) kk'K.' \coshp' cosh3p' cosh hu' 



(4.14) 



multiplication of ( 4.14 ) by k' gives (|4.11|) . 
We will discuss approximations 



£(°)( M ) = Asmhu'i (4.15) 
(u) = A sinh v! - Ai sinh 3m', (4. 16) 

£ {n) (u) = Asinhu' - A 1 smh3u + ... + (-l) n A n smh(2n + l)u (4.17) 



Within our notation A = Aq. We note, that u' and Ain+un — 0,1,2, have the 
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following expansions 



"' = ^" S3 ( 1 -i*° + -")" < 4 - 18 > 



^ = -{r Q k +- Q k +0 { k )) ( 4 - 2 °) 

A 3 = —k' 7 + O ( k' 9 ) (4.21) 
3 4096 V J K J 

However with our notations we do not expand Ai and u' in ( |4.18|) - ( |4.21| ) in 
power series of k' 2 . 

Approximations that are truncated series expansion of ( }4.11| ), namely approxima- 
tions ( |4.15| ), ( [4.16| ) are shown in Fig.^. B. A comparison of these approximations with 
ones obtained within the conventional perturbation theory are presented in Fig.|^.C 
and |2|.D. In Fig.^.B the first two approximations calculated according to the modi- 
fied perturbation theory are shown. It can be seen that the second approximation 
is an improvement of the first one. A comparison of the first approximations calcu- 
lated according to conventional and modified perturbation theories are presented in 
Fig.|2|.C. These approximations almost coincide with each other. However, the sec- 
ond approximation calculated according to the conventional perturbation theory and 
the second approximation calculated according to the modified perturbation theory 
differ significantly as it can be seen in Fig.^.D. The second approximation calcu- 
lated according to the modified perturbation theory is closer to the exact solution in 
comparison with the second approximation obtained according to the conventional 
perturbation theory. Next terms of the expansion are shown in Fig.^j. 



4.3 Bounds for elliptic cosine expansion 

Let us prove the following lemma. 

Lemma 1 Given fi > and < c < 1. Than the following inequality holds 



I cn(Qt — K,k) — £( n-1 ) (Qt) I < ~(k'] 



(2n+l)(l-c) 



n 



1,2, 



for all 



< / < — ^= In — , < k' < 1/2 



nV2 k 

here & n \u) stays for the approximation defined by 



(4.22) 
(4.23) 



Remark 1 Estimation ( \4-23j) is not trivial in comparison to the proposition that the 
series ( 4-11 ) converges with \%Xu\ < K and — 1 < k < 1. The reason for that is that 
the estimate according to allows to consider t — > oo with k' — > 0. 
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.4. 

Figure 3: (color online) The function of t cn(t — K, k) is shown with the red line with 
k = 0.999. The green line depicts a linear approximation. The orange, purple and coral 
dashed lines depict the first, the third and the fifth approximations calculated by formula 
(4.11) respectively dotted lines of cyan, aquamarine and ochre color depict approximations 
of the second, the fourth and the sixth orders calculated by (4.11). 



Proof is based on the estimation of the residual term 

7T ( sinh((2n + l)p'u") sinh((2n + 3)p'u") 



kK' { cosh((2n + l)p') cosh((2n + 3)p') 
by the geometric sequence sum 



7r e 



-(2n+l)p'(l-u") 



that is, 

Taking into account 



kK' 1 - e -2p'(W) 



<>~ 2p{1 - u } = exp(-2— exp 



2K' 



2K' 



the estimation ( f4.26| ) makes if 



cxp 



7T 



K 



2K' 



7 exp 



7T 



nt 



2K' 



< 1 



(4.24) 

(4.25) 
(4.26) 

(4.27) 
(4.28) 
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If this condition holds then 



d n < - I k' 



2n+l 



3 (2n+l)f2t 



< i^'^+iXi-e) ->■ 



A- 



(4.29) 

/or k' ^0 and c < 1 (4.30) 
The condition ( |4.28|) is satisfied if there is a restriction on t, such that 



V i 

k J exp (2f2t) < - 



(4.31) 



Taking into account the latter limitation and the property of the elliptic functions 
formulated below the condition ( (4.28| ) we get 



cxp 



7T 



K 



cxp 



7T 



2K' J 1 \2K' 
The detailed proof can be found in Appendix A. 



< k'e 



i„nt 



(4.32) 



5. n-mode approximation 

In the present section approximations of exact solutions are presented and accuracy 
of these approximations is estimated for small coupling constant e . 

Taking into account the series ( [4.1 ID let us consider an approximation of the 



exact solution (|3.18|) . We get 

oo 

q(t) = a cn(ttt - K, k) = ^ s n(t) 



n=0 



where the following notations are introduced 
'nttt 



So(t) = Asinh 



2K' 



S n (t) 



A = a 



71 



l) n A n sinh 

7T 



kK' , /ttK 
cosh 



2K' 



K = a kK' 



f(2n + l)7rfit 
1 



/(2n+l) 
COSh (^K^ 



7T 



K 



Let us introduce the notion of n-order mode approximations 

Qo(t) = s (t) 

qi(t) = s (t) + si(t) 

q 2 (t) = s (t) + si(t) + s 2 (t) 



(5.1) 



n = 1,... (5.2) 
, n = l, ...(5.3) 



(5.4) 
(5.5) 
(5.6) 



q n (t) = s {t) + si(t) + ... s n (t) 



(5.7) 
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We note that arguments of the sinh functions and A+, i = 1, ... are power series 
of k' 2 , however in definitions (|5.4j) - Q5.7T) we do make these expansions. 

In order to estimate an error of approximation of ( |5.4j ) - (|5.7| ) the following 
estimations for k' < 1/2 are used. The following Lemma provides the estimations 

Lemma 2 For k' < 1/2 the following estimations hold 

1. 

a = ^( 1 + T3Wr<7? <5 ' 8) 

2. 

, in 1 , 1 n 2cE 

^( 1+ ^) 1/4 <,( 1+ f) 

Proof (|5.8|) follows from the fact that if k' < 1/2, then (see. Fig. |4|. A) 

^r<2 (5.11) 

1 - 2k' 2 v ; 

Taking into account (|3.19|) we get (|5.^ ). 

The estimation ( f).9\) we get in the following way. The estimate holds, see Fig. ^ 

B, 

^- ( l -nrr=) <1 ^ o< y <oo (5.12) 




Figure 4: A. The estimate (5.12) is shown in the plot: k' 2 is along horizontal axis, and 
the left-hand side of (5.12) is on the vertical one. B. The illustration of the estimate (5.13). 
y is along the horizontal axis, and the left-hand side of (5.13). along the vertical one. C. 
The illustration of the estimate (5.15). y is along the horizontal axis, and the left-hand 
side of (5.15) is on the vertical axis. 



Therefore 

./o 1, 1 2eE 



I ^ l^-T- (5-13) 
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The estimate 

(1 + 1/) 1/4 
1+ y 

is shown in Fig.0 C. We get 



< 1 (5.14) 



< H 1+ ^ ) {r,Ar,) 



From the estimations ( |5.9| ), (|A.21 ) and ( |5.8| ) the estimation of the exact solution 
follows (O) 



\q(t)- qi (t)\ < adi < 2/^/-- 



14 /2e£\ 



3(l-c)/2 



Sv 7 ^ (2eE 



l-Sc 



for 



Let us note that according to (FT 



c 1 
0<t< ;= In — , c<l/2. 



//2 ^2 fc' 2 ' 



ona has 



1 > /i 4 



A;' 2 - 2e£ 



2e£ 



,1 i A* 
n — > In J — 

k 2 ~ 2eE 



i.e. if the following is true 

then ( p. 17| ) holds. 

In the same way we get 

Lemma 3 For 



c . a 
t < = In H 



^2y/2 2eE 



c , u 
t < t= In p 



estimations 



hold. 



\q{t) - q n -i{t)\ < 



^2^2 2eE 



, 2n + l 

8V2E (2eE\ " ~' 



kjj, \ jj, 



n = 1, 



Proof follows from (|A.25[) and estimations 



l<K*)-<k-i(f)l<y(*' 



(2n+l) 



e (2n+l)m < a ^/^,/\(2n+l)(l-c) 

k 



8V2E (2tE 



kfi \ /i 4 



(5.16) 

(5.17) 

(5.18) 
(5.19) 
(5.20) 



(5.21) 



(5.22) 



(5.23) 
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Explicit form of its lowest terms is 



n - 3 n - 75 / 



1 / =A 

16' 11 256' 



(6.11) 



Using approximation 

= (1 + fc ,2 Oi + ...A; /2 "O n ), (6.12) 
A (n) = A (k m l i0 + k' 2 ^l lX + ... + k'^ + % n ) . (6.13) 

we will obtain the following approximations for each term in the series ( |5 . 1| ) 

4 0) (t) = Asmh(fit), (6.14) 
s l o\t) = A sinh ((1 + k' 2 n x )nt) , (6.15) 
4 n) (^) = ^sinh ((1 + k' 2 nx + ...k' 2n n n )fit) . (6.16) 

We notice that, with adopted notations we do not expand A in power series of k' 2 . 
In the same way, 

s<™\t) = (-iy A\ n) sinh ((2i + l)fin®t) , (6.17) 

or explicitly 

s ^\t) = (-l)M (fc^o + ... + k' 2 ^% n ) sinh ((2i + + fc /2 fia + ....//'O,)/) 

(6.18) 

In accordance with ( |6.9| ) and ( 6.13 ) we have the following approximations for 
si(t) 

s?'°\t) = - A k'\ smh{3 (it), (6.19) 
si 0,1) (t) = -A k' 2 l lfi sinh (3/i(l + k' 2 n x )t) , (6.20) 
sf' 2) (t) = -A {k' 2 l lfi + fc /4 Zn) sinh (3/i(l + fc /2 fii + k' 4 n 2 )t) . (6.21) 

We will also use notations: 

sf xact ' j \t) = (-iy Ai sinh ((2z + + + ... + fc /2i %) t) . (6.22) 
s (n, e ,.a Ct)(t) = A (jfc/M^ + + k l2 ^% n ) sinh ((2z + l)/ifit) . (6.23) 

6.2 n-mode approximation with power series of A 

In order to express the series (|5.1|) in a way analogous to (|2.13|) we have to expand 
O and Aij A in power series at e 



f2 = ji 



A,=A 



1 + 



A 2 e^ 



Ah 



A 2 e 
T 2 " 



n 2 + 



t+1 



Oi = -, h = --. 

8 32 



(6.24) 

(6.25) 
(6.26) 
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We will also use notations 



(n) = fi 

A (n) = a 



1 + 



Ah^ 



Ah 
fj? 



k + 



Ah. 
» 2 



1* 



fl 2 + ... + 
hi + ...+ 



Ah 
Ah 



In the lowest terms we have the following approximations 



sf(t) 



A sinh ( fit) 
A sinh ((1 + 



Ah 



Qi) fit 



^(o,i) 



(t) 



-A 



A 3 e - ( A 2 e - 

l x sinh 3( 1 + — Oi) at 



(6.27) 

(6.28) 
(6.29) 

(6.30) 
(6.31) 

(6.32) 



Ah, MV 2 



h + 



sinh 3 1 



Ah^ (Ah^ 2 



-Oi + 



n 2 /it 



(6.33) 



here Ox and l\ are given by ( |6.26 ). The bar over Sj stays for expansions in power 
series of e (or A). 

We have the following Lemma 

Lemma 4 When amplitudes (\5.3j) and ( 2.19[) are equaled, i.e. 

21 = A (6.34) 

the following relations hold 



k 

rru 



Ink 

a 



Lemma (HI) allows to claim that if 21 = A then 



(t), 



ti j \t) = s^\t) n = l....,i,j = 0,l. 



(6.35) 
(6.36) 

(6.37) 
(6.38) 



7. Bounds to the hyperbolic analogue of the modified pertur- 
bation theory 

The goal of the discussions below is to obtain an approximation error estimate for 
( p,19|) mentioned at the beginning of the paper. To do this we mention that in 
approximation ( |5.7D constant Ai and the sinh functions argument are the power 
series of k' 2 , and consequently the powers of e. 
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7.1 One mode approximations 

Approximations to the exact one mode expression 



have the form 



</n(f) = .4sinh( (7.1) 



s ( o\t) = Asmh(fit) (7.2) 
$\t) = AsmhL(l + ^et)^ (7.3) 

S f) (t) = A sinh (V f 1 + % - ^eA t) (7.4) 



V 2 8 /z 



In relations ( |7.2| ) - ( [7.4| ) the sinh-functions arguments are different from ones in 
J]). We denote these differences by A n . (calculations in the file approx-ch.mw). 

A i = U U + ^ - — U (7.6) 



8/x 2 y 2K' 

Let us estimate contribution due to the mentioned differences. We have 

So(t) - 4 n) W = A [sinh (jiflt) - sinh (/iOt + A (n) )] (7.8) 
= A [sinh (/xf2 £) — sinh (/ifi t) cosh(A„) — cosh (/ifi t) ■ sinh(A n )] 
= A sinh (/xfl t) (I - cosh(A n )) - A cosh (fifl t) ■ sinh(A„) 

Taking into account that for small A, for instance, < A < 1, the following relation 
holds 

sinh(A) < 2A, (7.9) 
|l-cosh(A)| < 2A, (7.10) 

we get for t > and 

A n < 1, (7.1 1) 

\s (t) - sfr\t)\ < 2A n A sinh (fiQt) + 2A n A cosh (fiftt) 

= 2A n Ae fJ,nt , (7.12) 

(7.13) 
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At this point one can see that the error of the approximation of Sq (t) to So(t) 
is defined by infinitesimality of A n . 
In particularly, for Aq we have 



A = fitSo, 



5n = 1 



7T 1 



2K' VI - 2k' 2 



From analysis (see Fig. |5|) we get 



and therefore, 



S < V2k' 2 , < k' < -, 



A < V2fik'H, < k' < ]- 



(7.14) 
(7.15) 



(7.16) 



(7.17) 



0,3 0,4 0,5 




Figure 5: The estimate (7.15) is showed. Along the horizontal axis k' 2 are plotted, and 
on vertical axis the relation the relation of the right-hand side of (7.15) to k' 2 . 



Estimates A n , n — 1, 2 are presented in Appendix. B. 

7.2 Two modes approximations 

Consider the two modes term qi(t), 

qi{t) = q (t) + Sl (t), 

where 



si(t) = — Ai sinh 



A\ = a 



7T 



3nQt 

2K' 
1 



n = 1, 



kK 1 , /3ttK 

cosh 



2K' 



(7.18) 

(7.19) 
(7.20) 
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and as an approximation to Si(i) take 



S i°' 0) = -^sinh(3/it) (7.21) 

or 

sf'°\t) = -A k' 2 ^- sinh(3//t) (7.22) 
16 

A deviation of from the exact si(t) is not only due to changes in the argument 
of sinh - function, 3Ao, but also due to corrections to the exact value of A\, 

A 1 = ^k' 2 A + a k' 5 L 1 (k r ) (7.23) 

and we have an estimation 

Lt < 0.1 for k' < 1/2. (7.24) 
This estimation is true since (see Fig. |6]). 

k' 2 IT 1 7T 1 1 5 

"TefcK 7 " /ttK\ + fcK 7 " /3vrK\ < 10 " (7 ' 25) 
cosh — — cosh 



2K' 



/3ttK\ 

A 

V 2K ' / 



We get 



sS 0,0) (t) - sUt) = -A k' 2 — sinh(3ut) + A x sinh(3ufit) (7.26) 

16 

= -A k' 2 — sinh(3/i(5 Q + + A x sinh(3/ifit) (7.27) 
16 

= -A k' 2 — (sinh(3/i(<5 + Q)t) - sinh(3//ft£)) (7.28) 
16 

k' 2 

+ o( A + Ai) sinh(3/ifit) (7.29) 

16 



Consequently 



\sf' 0) (t) - si(t)\ < -Ak' 2 — | sinh(3/i(5 - «)*) - sinh(3//fit)| (7.30) 

16 

+ a— | sinh(3/xfi)t| (7.31) 



Taking into account an estimation of the form (|7.12 ) we get, 

|sf' 0) (t) ~si(t)\< Ck' 4 e 3 ^' (7.32) 

where C is constant. 
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Figure 6: (color online) A. The estimation (7.25) is illustrated. On the horizontal axis 
k' 2 is plotted, on the vertical axis the values of the left-hand side of (7.25) (red color) and 
the right-hand side (7.25) (green color) are presented. B. The zoom of the plot in the left 
panel A for small values of k' 2 . 



7.3 The second order of the modified perturbation theory approximation 

The following theorem holds 
Theorem 1 For 

0<t<— ^ln-^— and 0<c<- (7.33) 
2V2fM 2eE 2 y 1 

there are C and a < 3/4 such that, the following estimation is true 

m -4»(t)- s r\<C?^( e 4) 2 ~°, „=!,... (7.34) 

Remark. In the inequality (17.34) C - dimensionless constant, factor is 

H 

introduced to get dimensionless quantity, - dimensionless factor. 

m 

Proof. It follows from the estimates (|5.22|) and ( [7.32|) . Indeed, from (|5.22 ) it 
follows that 

2-#c 



8V2£ ( 2eE \ ' -' 

-Qi{t)\< 

In accordance with (15.51 



\q(t) - gi (t)\ < ^— i^— ) . „ = !.... (7.3:)) 



Qi{t) = q (t) + si(t) (7.36) 

and in the estimation (|7.34j) we took the approximation s^\t) to qo(t) and the 
approximation s[°'°' ) to Si(t). We get 

q(t) - s$\t) - 4°' 0) = q(t) - qi(t) + q (t) - s$\t) + 8l {t) - s?>°\t), (7.37) 

consequently 

\q(t) - s$\t) - 4°' 0) | < \q(t) - qi (t)\ + \q (t) - s£\t)\ + \ Sl (t) - 4°' 0) (t)| (7.38) 
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In order to finish the proof we need the estimates ([7.11]) and ( B.9 ) for 
Gathering these estimations we get ( |7.34[) . 



7.4 Bound of the n-mode approximation 

Consider an n-mode s n (t), 

f(2n + l)irttt 



s n (t) = (-l)M n sinh 



V 2K' 

l)MZ n sinh((2ra + l)/zfit), n=l,... (7.39) 



where 



* 1 

A n = a 



kK' , f(2n + 1)ttK 

cosh 

V 2K' 

/vrK 

cosh 



V 2K'7 Acosh(KfJ) 



/(2ra+l)7rK\ cosh((2n + l)Kfl) 
C ° S ^ 2~K' ) 



and 



n[ ' cosh((2n + l)Kft) 1 } 
As an approximation to s n (t) we take 

s (r,m) = (_^n A j(r) ginh ^ + ^qHj) ( 7>42 ) 

where the series Zn ' is obtained from the expansion of the right-hand side of ( |7.41|) 

1$ = ln, k' 2n + ...+ln,rk' 2{n+r) (7.43) 

and O 1 " 1 ) is defined in (|B~9D, 

ft (m) = 1 + n x k' 2 + ... + O m fc ,2m (7.44) 
The following lemmas are true . 



Lemma 5 For < k' < 1/2, r > l,n > 1, t/iere are constants and L^h 1 such 
that 

in(jfe') - (r) (/c')| < B {r) k' 2{r+1) (7.45) 

|Z n (A/) - < 4 r) ^ /2(n+r+1) . (7.46) 
From Lemma [5] one gets the following 
Lemma 6 



n (t) - 4 r ' m) WI = S^k' 2 ( n+r+1 >e^' m)t (7.47) 
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7.5 Error of the n-mode approximation by power series of A 

The n-mode can be considered function of A: 

s n (t) = (-l)M Z n (A) sinh ((2n + l>n(A)t) , n = 1, ... (7.48) 
Using the power series on A one introduces notations 

g(r,m) = A J(r) ^ ^ + ]_)^fjM^ (749) 

where and f2( m ) are obtained via reexpansions of the right-hand side of ( [7.41] ) 
and ( [7.440 in a power series of A 

W = l nfi \ n + ... + l n , r \( n+r \ (7.50) 

fV m > = l + n 1 A 2 + ... + O m A Tn (7.51) 
The following lemmas are true. 

Lemma 7 For A small enough there are constant and Ln such that 

|0(A) - O w (A)| < fiWA (r+1) (7.52) 

|S»(A)-i£>(A)| < 4 r) A (n+r+1) - (7.53) 
From Lemma || it follows that 
Lemma 8 

|s n (t) - 4 r ' m) WI = ^ r ' m) A (n+r+1) e^ r ' m)t (7.54) 



The following theorem holds ||18 
Theorem 2 For 

< t < 



, A* , 1 

ln — — and < c < - 

2e£ 2 



one can /md constants C n and o < 3n/4 such that, the following estimate is true 

2n—a n 



(7.55) 



q{t)-q {n \t)\<C ri 



E ( eE 



n 



1. 



where 



q {n \t) = j2sr t ' n ~ v (t) 



-(n—i,n—i) , 



(7.56) 



(7.57) 



i=Q 



In £/ie right-hand side of (\7.57 ) contributions to all the first n-modes are taken into 
account, and in each of them approximations to the exact values of the amplitudes 
and the phase are taken with orders not exceeding (n — i) 
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v E 

Remark. C n is a dimensionless constant, the factor is introduced by the 

dimension analysis, ^ is a dimensionless factor, a n depends on n. 

Proof follows from the estimates ( 5.22 ) and (|7.32|) . Indeed, from ( 5.22 ) it follows 

that 



2n+l . 



\q{t)-q n {t)\ < 



8V2E f2tE\ 2 2 



n = 1 



(7.58) 



According to 

q n (t) = s (t) + Sl (t) + ... + s n (t), (7.59) 

and in ( |7.57|) we take the approximation of sjj (i) to qo(t) and the approximation 
1,n ~ , for si(t) and so on. We get 

q(t) ~ q (n) (t) = q(t) - s (t) 



■ ■ - s n (t) (7.60) 

+ 8 (t) ~ t\t) + Si(t) - st'^it) + .... + S n (t) - 4°'°)(t), 



and consequently 



\q(t) - qW(t)\ < \q(t) - q n (t)\ + \q (t) - $>(t)\ 

+ \ Sl (t) - S f - 1 '- 1) (t)| + ... + \s n (t) - 4°' 0) (t)| (7.61) 

In order to finish the proof we need estimates ( [7.1 1|) and the estimate ( [B.9p for 
Having taken into account these relations we get ( |7.56| ). 



8. Conclusion 

The method for approximate solution to nonlinear dynamics equations in the rolling 
regime is presented. It is shown that in order to improve perturbation theory in the 
rolling regime it turns out to be effective not to use an expansion in trigonometric 
functions as it is done in case of small oscillations but use expansions in hyperbolic 
functions instead. In particular the Higgs equation in the rolling regime is considered. 
This regime is investigated using the representation of the solution in terms of elliptic 
functions. An accuracy of the corresponding approximation is estimated. 

In this paper we have investigated the rolling regime for motions starting from 
the top of the Higgs potential, see Fig.|I].A. The same method can be used for motion 
starting with zero velocity from the side that is near to the top, see Fig]l|.C. 

As to possible cosmological applications of rolling solutions it is worth to mention 
that in cosmology various rolling solutions are widely used. There is a notion of the 
slow roll regime that means that one can ignore the second order derivatives in the 
equation of motion J7|, The rolling from the top of the potential is not necessary 
described by the slow roll approximation and methods similar to ones developed in 
this paper seems to be rather suitable [[HJ . 



-27- 



There are rolling solutions also in the nonlocal cosmology [19. 20, 21, 22, 23 
These solutions have their analogue in the flat space-time [25|, ^6]. Non slow 
roll regime is important in inflation and it is related with so-called stretch effect in 
nonlocal theories [E7, I 
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A. Proof of Lemma 1 

Consider first in the case n — 1. The difference in the left-hand side of ( |4.22| ) is 

sinh 
cosh 



A, = cn(nt — K,k) — JL ( A. I ) 



2K' 



Using (14.11 ) and the majorant of the following sum 



7T f sinh(3 P V') sinh(5pV') 1 

1 jfeK' 1 cosh(3p') cosh(5p') "J 1 ' J 

where 

,, fit , 7rK ,, , , nQt 



we get 

i r I < 71 V" sinh((2j + l)p'u") 7T e _ p / ( i- u ») 2 . 7 y(i-u") 
1 itK'^ cosh((2j + iy) jfeK' ^ 

Using 

fcK' ^ " jfeK' 1 - e-V(W) ' 

we get 

|n| < d t , (A.5) 



where 



7T 6 



In Fig.0. A. the bound ( |A.5| ) is shown. 



- 28 - 



We notice that the series in (|A.4|) converges if 

e -p'(W) < 1; 

that is true if 

< u" < 1, 



since 



7rK 

2K 7 



> 



nt 



Since u" = — the relation (|A.8 ) means that 



Qt 

0<-<l 



Consider the right-hand side of (|A.5| ) 

n e -3p'(l-u>) 



di 



71 



7rK^ 


exp 




2K'; 


V 2K' J 


7rK\ 

- 2 2K' 


exp 


y, 2K' 



(A.7) 
(A.8) 

(A.9) 



(A.10) 



(ah; 






4 A. 



Figure 7: A. The estimate (A. 5) is shown. Along the horizontal axis time values t are 
plotted (with = 1), and on vertical axis the relation ^ is plotted. B. The estimate 
(A. 12) is shown. Along the horizontal axis time values k are plotted, and on vertical axis 
L(k) given by the relation (A. 12) is plotted. C. The estimate (A. 13) is shown. Along the 

7T 

horizontal axis time values k are plotted, and on vertical axis is plotted. 
Using explicit estimates (see. Fig. [7].B and Fig.[7].C, respectively) 

7rK 



L(k) = exp 



2K' 



y/l - k 2 < 



7T 



2K' 



< 1 



we have 



exp 



7T 



K 



2K' 



exp 



2K 7 



< fc'e 



(A.12) 
(A.13) 

(A.14) 
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We require that 



to estimate flA.ll ). 



Then we have 



kK 

Notice that for kl < 1/2 



^k'^j exp(2Qt) < ^ (A. 15) 



dl < 2^- (k f ) \ 3nt < \ fk') 3 e 3nt (A.16) 



and ( |A.15| ) are true if 

(k'^j exp(2Qt) < (k'^j exp (2V2 fitj < ^, k' exp (y/2fjti) < -j= (A.18) 



Consequently (|A.18|) is true, for 



*< -^Inl/Jf, c<l/2. (A.19) 



Indeed. 



exp (v^/xt) < (A;') 1_c < l/v 7 ^ if c < 1/2, k' < 1/2. (A.20) 
With ( |A.19| ) the following is true 



tZi<iM 3 e 3nt <^(ife , ) 3(1 - c) . (A.21) 



In the same way for 



vr f sinh((2rc + l)pV') sinh((2n + 3)pVQ } 
r ™ jfeK' \ cosh((2n + l)p') cosh((2n + 3)p') J 1 ' ' 

we have 

|r„| < dn, (A.23) 

where 

^ e -(2n+l)p'(W) 



and 



rfn " jfeK' 1 - e-V(i-«") ' (A ' 24) 



A f \ 2n+1 A 
d n <-(k') e ^n+mt < 4 (A ,/ )( 2n+l)(l-c)_ (A>25) 

k\ I k 



For small fc' < 1 /2 we get k > 1/2 and 
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(li <-ik' \ e 3nt <^(k') e 3nt 



Finally for k' < 1/2 k > 1/2 we get 



4 



di<~{k') e 3nt <8(k') e 3nt 



B. Bounds of A 



Having introduced notations 
we obtain: 



5 = 1 

51 = 



71 



2K' VI - 2k' 2 

3A 2 \ 7T 

1 + — — e 



1 + - 1 + 



2K' VI - 2k' 2 
{ 



1 - 2k' 2 
3A 2 15A 4 o 



7T 

fcK 7 



7T 



K 



1 + 



v cosh V2ic;/ 

7T 1 



7T 



2K' VI - 2k' 2 



8/i 2 2 V 2K' VI - 2k' 2 



( 



i + ; i 



1 - 2A;' 2 



fcK' 



v cosh (ii;/ 



15 

2« 



1 - 2A;' 2 



7T 



1 



2K' VI - 2A;' 2 
These formulae are based on the relations: 

TlQ 7T [I 



2K' 2K' VI - 2A;' 2 



/ 



/' 2 



We have (see Fig. ||) 



1 



1 - 2k' 2 



\ 



fcK' 



7rK 



v cosh \2^;/ 



8 < V2k' 2 , < k' < - 



(A.26) 



(A.27) 



(B.l) 

(B.2) 
(B.3) 



(B.4) 

V 



fcK' 



7rK 



v cosh V2K';/ 



(B.5) 



(B.6) 



(B.7) 



31 



and, consequently, 



We also have 



A <V2fik' 2 t ,0 <£'<-. 



8 X < k'\ 

This estimation is illustrated in Fig.|[A. 



(B.8) 
(B.9) 





0,2 0.3 



Figure 8: (color online) A. The plot illustrates the estimation (B.9). On the horizontal 
axis k' 2 is plotted, on the vertical axis the values of the left-hand side of (B.3) (red color) 
and the function A:' 4 (green color) are presented. B. The plot illustrates the estimation 
(B.4). On the horizontal axis k' 2 is plotted, on the vertical axis the values of the left-hand 
side of (B.4) (red color) and the function k' e (green color) are presented. 



From ( p.9|) we get that 

Ai < fik'H. (B.10) 

In Fig||B. we show functions 62 = ^(fe') and k' 6 . From from this plot one can 
see that 

5 2 < k l& (B.ll) 

and 



A 2 < fik'H. 



(B.12) 



- 32 - 



References 



[1] N. Krylov, N. N. Bogoliubov, Introduction to Nonlinear Mechanics. Princeton 
University Press, Princeton, N. J., 1947. 

[2] N.N. Bogolyubov, Y. A. Mitropolski, Asymptotic Methods in the Theory of 
Nonlinear Oscillations. Gordon and Breach, New York, 1961; 

[3] V.I. Arnold, V.V. Kozlov, A.I. Neishtadt, Mathematical aspects of classical and 
celestial mechanics, Springer, 2006. 

[4] V. V. Kozlov, S. D. Furta, Asymptotic expansions of solutions of strongly nonlinear 
systems of differential equations, Izhevsk, Regular and Chaotic Dynamics, 2009. 

[5] V. A. Rubakov, Classical theory of gauge fields, Princeton University Press, 2002 

[6] A. Linde, Inflation and Quantum cosmology, Academic Press, Boston, 1990. 

[7] V.F. Mukhanov, Physical Foundations of Cosmology, Cambridge University Press, 
2005. 

[8] D.S. Gorbunov and V.A. Rubakov, Introduction to the Theory of the Early Universe, 
URSS, Moscow, 2009 (in Russian). 

[9] I. Y. Aref'eva and I. V. Volovich, "Cosmological Daemon," JHEP 1108, 102 (2011), 
arXiv:1103.0273. 

[10] I. Y. Aref'eva, N. V. Bulatov and R. V. Gorbachev, "FRW Cosmology with 
Non-positively Defined Higgs Potentials," arXiv:1112.5951. 

[11] IV. Volovich, Randomness in classical mechanics and quantum mechanics, Found. 
Phys. 41:3 (2011), 516-526. 

[12] IV. Volovich, Bogoliubov equations and functional mechanics, Teoret. Mat. Fis. 
164:3 (2010) 354-362. 

[13] EV. Piskovskiy, IV. Volovich, On the Correspondence between Newtonian and 
Functional Mechanics, Quantum Bio-Informatics IV, eds. L. Accardi, W. 
Freudenberg, and M.Ohya, World Sci. (2011) pp. 363-372. 

[14] A.M. Jouravsky, Handbook of elliptic functions (in Russian), AS USSR, 
Moscow-Leningrad, 1941 (in Russian). 

[15] M. Abramowitz, and LA. Stegun, Handbook of mathematical functions: with 
formulas, graphs, and mathematical tables, Dover Publications, 1964 

[16] N.I. Akhiezer, Elements of the theory of elliptic functions, AMS, Providence, RI, 
1990. 

[17] IYa. Arefeva, EV. Piskovskiy, and IV. Volovich, Rolling in the Higgs Model and 
Elliptic Functions, submitted to Theor. Math. Phys. 



- 33 - 



[18] I.Ya. Arefeva and I.V. Volovich, Asymptotic expansion of solutions to the rolling 
problem, submitted to Proc. of Steklov Mathematical Institute. 

[19] I.Ya. Arefeva, Nonlocal String Tachyon as a Model for Cosmological Dark Energy, 
2006 AIP Conf. Proc. 826, pp. 301-311, astro-ph/0410443 

[20] I.Ya. Arefeva and L.V. Joukovskaya, 2005 Time lumps in nonlocal stringy models 
and cosmological applications, JHEP 0510, 087, hep-th/0504200 

[21] G. Calcagni, Cosmological tachyon from cubic string field theory, JHEP 0605, 012 
(2006), hep-th/0512259 

[22] I.Ya. Arefeva, A.S. Koshelev, Cosmic acceleration and crossing of w = — 1 barrier in 
non-local Cubic Superstring Field Theory model, JHEP 0702:041, 2007, 
hep-th/0605085 

[23] I.Ya. Arefeva, L.V. Joukovskaya, and SYu. Vernov, Bouncing and accelerating 
solutions in nonlocal stringy models JHEP 0707 (2007) 087, hep-th/0701184 

[24] V.S. Vladimirov and Ya.I. Volovich, "On the nonlinear dynamical equation in the 
p-adic string theory," Theor. Math. Phys. 138 (2004) 297, math-ph/0306018. 

[25] Ya. Volovich, "Numerical study of nonlinear equations with infinite number of 
derivatives," J. Phys. A A 36, 8685 (2003) math-ph/0301028 

[26] I.Ya. Arefeva, L.V. Joukovskaya, and A.S. Koshelev, Time evolution in superstring 
field theory on nonBPS brane. 1. Rolling tachyon and energy momentum 
conservation, JHEP 0309 (2003) 012 

[27] N. Barnaby, T. Biswas and J.M. Cline, "p-adic Inflation," JHEP 0704 (2007) 056 
hep-th/0612230 

[28] J.E. Lidsey, Stretching the inflaton potential with kinetic energy, Phys. Rev. D76, 
043511 (2007), hep-th/0703007. 

[29] N. Barnaby and J.M. Cline, "Large Nongaussianity from Nonlocal Inflation," JCAP 
0707 (2007) 017, arXiv:0704.3426 

[30] G. Calcagni, M. Montobbio and G. Nardelli, "Route to nonlocal cosmology," Phys. 
Rev. D 76, 126001 (2007), arXiv:0705.3043 

[31] L. Joukovskaya, Dynamics in Nonlocal Cosmological Models Derived from String 
Field Theory, Phys. Rev. D 76, 105007 (2007), arXiv:0707.1545 

[32] I.Ya. Arefeva, L.V. Joukovskaya, and SYu. Vernov, Dynamics in nonlocal linear 
models in the Friedmann-Robertson-Walker metric, J. Phys. A: Math. Theor. 41 
(2008) 304003, arXiv:0711.1364 

[33] N. Barnaby and J.M. Cline, "Predictions for Nongaussianity from Nonlocal 
Inflation," JCAP 0806 (2008) 030 arXiv:0802.3218 



- 34 - 



[34] L.V. Joukovskaya, "Dynamics with Infinitely Many Time Derivatives in 

Friedmann-Robertson- Walker Background and Rolling Tachyon," JHEP 0902 (2009) 
045, arXiv:0807.2065 

[35] D.J. Mulryne and N.J. Nunes, Non-linear non-local Cosmology, AIP Conf. Proc. 
1115 (2009) 329-334, arXiv:0810.5471 

[36] N. Barnaby, Nonlocal Inflation, Can. J. Phys. 87 (2009) 189-194, arXiv:081 1.0814 



- 35 - 



